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Abstract 

Methods to estimate the Hausdorff dimension of invariant sets of 
scattering systems are presented. Based on the levels' hierarchical 
structure of the time delay function, these techniques can be used 
in systems whose future-invariant-set codimensions are approximately 
equal to or greater than one. The discussion is illustrated by a nu- 
merical example of a scatterer built with four hard spheres located at 
the vertices of a regular tetrahedron. 
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Usually, in a scattering process almost all particles outcome asymptoti- 
cally free after interacting with a potential. Some particles, however, remain 
in the scattering region bouncing forever under the action of the potential. 
The nature of these trapped orbits is in general very intricate, including 
fractal dimensions and chaotic behaviour. The dimension seems to be the 
most fundamental quantity characterising this set. When the dimension is 
large enough, scattered particles may present sensitive dependence on ini- 
tial conditions reminiscent of the chaotic dynamics of the trapped orbits. 
This phenomenon allows the computation of the uncertainty dimension of 
the bounded orbits, defined as follows. For a scattering system, colour by 
black initial conditions of particles scattered, for instance, to the right-hand 
side of the scattering centre and by white initial conditions corresponding 
to particles scattered to the left-hand side. An initial condition of a given 
colour at a distance less than e of initial conditions of a different colour is 
labelled e-uncertain. The final state (left or right) of a particle associated 
to an e-uncertain initial condition is uncertain if this point is determined 
only within a precision of order e. In the limit of small e, the fraction of 
such ^-uncertain points is f{e) ~ e^~^^, for denoting the dimension of 
the subspace where the initial conditions are taken, and Du the correspond- 
ing uncertainty dimension of the intersection of the trapped orbits with this 
subspace. Therefore can be computed by fitting In /(e) as a function 
of ln(e). This procedure works in systems where the set of bounded orbits 
has dimension greater than — 1 in the A^-dimensional subspace of initial 
conditions. That is supported by the fact that in this case is expected to 
be equal to a more general concept of dimension, the Hausdorff dimension. 

The Hausdorff dimension of a set A is defined as follows Let Sj 
denote a countable collection of subsets of the Euclidean space such that the 
diameters ej of the Sj are all less than or equal to 6, and such that the Sj 
are a covering of A. We define the Hausdorff sum 



where the infimum is taken on the set of all the possible countable coverings 
of A satisfying ej < 6. Then we define the d- dimensional Hausdorff measure 



(1) 



K{d) = \imKs{d). 







(2) 
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It can be shown that K{d) is infinite if d is less than some critical value 
Dh, and is zero if d is greater than Dh- This critical value is the Hausdorff 
dimension of the set A. 

Among the several existing concepts of dimension, the Hausdorff dimen- 
sion is the one that better implement the intuitive notion of dimension that 
we have in mind. Unfortunately, the uncertainty methods cannot be used to 
estimate Dh when its value is less than or equal to — 1. Even when Dh 
is slightly greater than — 1, its evaluation through the numerical compu- 
tation of the uncertainty dimension is highly inaccurate. The aim of these 
work is to present methods to compute the Hausdorff dimension of the set of 
trapped orbits - the future invariant set (//) - when uncertainty methods do 
not apply. The determination of this dimension is important, for instance, to 
understand the onset of chaos in scattering systems when a system parameter 
is varied. 

First in this communication, we introduce a method to take into ac- 
count the infimum requirement appearing in the definition of the Hausdorff 
dimension, that allows the direct computation of Dn^If)- Next, the Haus- 
dorff dimension is estimated from the box-counting dimension computed by 
using a proper-interior-maximum- /2A;e procedure. Finally, we define a new 
dimension, the self- similarity dimension, which is also used to approximate 
Dnilf)- Our methods can be applied for Dnilf) > A^ — 1 as well as for 
Dnilf) < N — 1. The discussion is illustrated along the text by a numeri- 
cal example of a simple scattering system presenting Dnilf) < N — 1 and 
Dnilf) > N — 1 for two ranges of values of its control parameter. 

The example we shall consider is due to Chen et al. , who investigated a 
system consisting of a particle bouncing between four identical hard spheres 
placed at the vertices of a regular tetrahedron of unit edge. The centres 
of the spheres are located at the coordinates (0, 0, v^/3), (1/2, — 1/2\/3, 
0), (—1/2, — 1/2-\/3, 0) and (0, l/-\/3, 0), and their radius R is the variable 
parameter of the system. These authors contributed with a procedure to 
directly evaluate the Hausdorff dimension Dh{R) on 1-dimensional lines. By 
taking lines on the plane P = {z = —k {k > R),Px = Py = 0,pz = 1}, they 
computed the Hausdorff dimension Dh{R) of the intersection of // with P, 
obtaining Dh^OAS) = 1 + D//(0.48) = 1.4. It is reasonable to presume that 
Dh is an increasing function of R. Here we are interested in computing Dh 
for small values of the radius, where the dimension of If (1 P is expected to 
be less than one (codimension greater than one). 
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We consider initial conditions {xq, yo) on P, and we define the time delay 
function T(xo, yo) as the number of collisions of the particle with the spheres. 
We also introduce C„ = {{xo,yo) \ T{xo,yo) > n}. Ci is the union of 4 discs, 
the projection on P of the four spheres. C2 is the union of 4-3 = 12 deformed 
discs, the projection on P of the reflection on each sphere of the other three 
spheres. C3 is the union of 4 ■ 3^ = 36 deformed discs, and so on. Therefore, 
for sufficiently small R, Cn is the union of 4 ■ 3"~^ deformed discs. The 
construction of the Cantor structure of the set // fl P = flJ^^Cn can be 
followed in Fig. 1, where we show Ci, C2, C3 and C4 for R = 0.37. 

Labelling the spheres by 1, 2, 3 and 4, the left shift symbolic dynamics of 
the future invariant set consists of sequences S1S2S3... of 1, 2, 3 and 4, where 
Sn is the number of the sphere of the n*'^ collision. The particle cannot collide 
two successive times with the same sphere, and hence the only constraint is 
Sn+i 7^ Sn (at least for R small enough). A relevant point here is that each 
one of the 4 ■ 3"~^ deformed discs of C„ can be identified by a sequence of 
collisions, c„j = [siS2...Sn] for j = 1,2, ...,4 ■ 3"~^. It allows us to define the 
Hausdorff sum 

Knid)= Y: |Cn,,f, (3) 

where \cn,j\ is the diameter of c„j, |c„j| = sup^^^g^^^ \x — y\, which can be 
easily estimated from a sample of initial conditions. The (i— dimensional 
Hausdorff measure K{d) is supposed to be the hmit of Kn{d) for n ^ 00. 
K{d) is infinite for d less than the Hausdorff dimension D^, and is zero for 
d greater than [^. Since it is not known how to compute lim„^oo-^n('^) 
numerically, we generalise an argument due to For n sufficiently large, 
the sums Kn{d) for different values of n will all intersect with each other at 
approximately the same point d = Dh- It provides a method to estimate the 
Hausdorff dimension. 

In Fig. 2 we show In K^id), \nK(i{d), \nKr{d) and InKgid) for R = 
0.37. The curves in fact intersect approximately at the same point giving 
Dh{0.37) = 0.87 ±0.01. In Fig. 3 (boxes) we show a graph of Dh computed 
by this method for 0.37 < R < 0.48 (it includes regions with Dh > 1 and 
with Dh < 1)- In the same figure (stars) we also show the box-counting 
dimension D^ of fl^^j^Cm. Dc was calculated by approximating nj^^^Cm by 
Cn for a sufficiently large n. The number of squares N{e) needed to cover C„ 
was counted for different values of the edge length e, and the dimension Dc 
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was then obtained by fitting InA^(e) = —Dclne. The comparison between 
Dc and Dh in Fig. 3 suggests that Dc{R) = Dh{R)- 

Some remarks about the computation of the box-counting dimension have 
to be made. The accurate evaluation of Dc requires the use of relatively 
small values for e. On the other hand, C„ has to reflect the fractal property 
of n^^^Cm, condition satisfied if |c„j| < e for all j. It implies the need of 
large values for n (we used n = 15 in our computations), demanding great 
computational efforts in finding C„ when is less than one. The remedy is 
based on the use of a modified PIM (proper-interior-maximum) procedure^ 
[^]. In the first step, we start with a sample of initial conditions on P taken 
at the corners of a grid. From this set we store those points that are in Ci as 
well as their first neighbours. We denote the stored set by Ci. In the second 
step, we reduce the edge length of the grid and iterate only those points of 
the new grid that are in the region covered by Ci. Then we store the points 
that are in C2 together with their first neighbours. We denote this stored set 
by C21 and so on. After n steps we obtain the desired approximation for Cn- 

Back to Fig. 1, in the centre of the figure we see an approximately self- 
similar structure. Each step of the construction of the Cantor set divides one 
disc in three smaller deformed discs. If we assume that the self-similarity is 
exact, we can introduce Ai, A2 and A3, the hnear scales by which the discs 
are reduced in one step. Under such hypothesis we have that Hausdorff and 
box-counting dimensions are equal and satisfy ||^ 

Af + Af + Af = l. (4) 

Even through the self-similarity is not exact in our case, (^) is expected to 
provide a good estimate for Dh (and Dc)- Fig. 1 also shows that the factors 
by which the discs are reduced are approximately the same. It allows us to 
take Ai = A2 = A3 = A, and we refer to the resulting dimension 

as the self-similarity dimension. The computation of A can be performed 
by a Monte Carlo method. We determine the areas of C„ and Cn+i, for a 
sufficiently large n, by evolving a sample of random initial conditions. An 



This technique can also be employed in the computation of |cnj|, to calculate Dh- 
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estimate of A is then given by 



An 



1 area(C^+i) 
^3 area(C„) ' 



A better statistics is obtained by taking A as the average on A„, A„+i, \n+k 
for an adequate k. For the tetrahedron scatterer, the value of A is very stable 
in relation to n and k, what allows us to compute a well defined A with 
not much computational effort. The resulting dimension is shown in Fig. 3 
(triangles) for 0.36 < R < 0.49 and is consistent with Dh and Dc computed 
directly from their definitions. We believe that Dg can be used to estimate the 
common value of Dh and Dc in most systems where some kind of statistical 
self-similarity takes place. A rigorous statement in this direction requires, 
however, further investigations. 

The dimension (D) that we computed corresponds to the intersection of 
the stable manifold of the invariant set (repellor) with the plane P. In the 
5-dimensional energy surface, the dimension of the stable manifold will be 
ds — S+D. The time reversibility of the dynamics implies that the dimension 
of the unstable manifold du will be equal to ds- Since the invariant set is 
the intersection of the stable and unstable manifolds, its dimension will be 
di = ds + du — 5 = 2D + 1. Therefore, the dimensions dg, du and di can 
be calculated by computing D with the techniques employed here. In this 
procedure, the critical value of the radius (Re) where D — 1 [Re ~ 0.41, see 
Fig. 3) does not play any particular role. In this context, to claim that the 
dynamics of the scattering systems is chaotic for D > 1 (rfj > 3) and regular 
for D <1 {di <?)) sounds largely a matter of semantics. What is sure is that 
there is a chaotic set of trapped orbits whose dimension increases with the 
radius R. The presence of this chaotic set prevents the existence of integrals 
of motion besides the energy (the system is non-integrable). 

In the computation of Dh-, D^ and Dg, we explored the hierarchical struc- 
ture of levels of the time delay function in order to infer the dimension of 
its set of singularities (flJ^^j^C^). The dimension of the singularities of the 
time delay seems to be the most fundamental quantity that we can measure 
to characterise invariant sets of scattering systems. The result involves no 
ambiguity since it will lead unequivocally to the dimension of the future in- 
variant set. The same is not true for basin boundaries, for instance. See Fig. 
4, where we show the basins of the tetrahedron scatterer defined by particles 
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scattered to x — +00 and particles scattered to x — >• —00, for R = 0.48 
(a) and R = 0.37 (b). The basin boundaries present smooth 1- dimensional 
parts on arbitrarily fine scales. These smooth parts are not on the future 
invariant set and, although immaterial when R = 0.48, they do affect the 
basin boundary dimension when R = 0.37. In this case the dimension of the 
fractal part {If fl P) is less than the dimension of the smooth parts. Smooth 
parts also forbid the use of uncertainty methods to estimate the dimension 
of the future invariant set from the time delay function. In computing the 
uncertainty dimension of the singularities, the uncertainty method effectively 
computes the box-counting dimension of the set of all discontinuities. The 
time delay is discontinuous on the (1-dimensional) frontiers of all C„ (see 
Fig. 1), what implies that the result will be always greater than or equal to 
1. The same happens with scattering functions. 

Finally, we remark that methods to directly estimate the Hausdorff di- 
mension of attractor sets (taking into account the infimum in (|I|)) was already 
discussed (see p and references therein). The main difficulty in adapting 
these methods to the case of nonattracting sets is that they require a sample 
of points on the set we want to measure. Even though a simple matter in the 
case of attractors, the computation of these points is nontrivial for nonat- 
tracting sets. A possibility would be the use of the PIM-triple procedure 
[^] to find trajectories which stay near the invariant set for arbitrarily long 
periods of time. The resulting methods are, however, far more complicated 
than those introduced here. 

In summary, we have presented, through an example, methods to esti- 
mate the Hausdorff dimension of nonattracting invariant sets of scattering 
systems. The methods apply, albeit not only, to systems with invariant set 
dimension so smaller that the uncertainty methods do not work. It includes 
cases where the codimension of the future invariant set is approximately equal 
to or greater than one. We stress that these methods are general and can be 
used in phase spaces of any dimension. And since they allow, in principle, 
the computation of arbitrarily small dimensions, these techniques should be 
useful to study how chaotic scattering comes about as a system parameter is 
varied . In particular, it may be interesting in investigations about routes 
to chaos in three dimensional scattering. 

The authors thank Fapesp and CNPq for financial support. 



7 



References 

[1] E. Ott, Chaos in Dynamical Systems (Cambridge University Press, New- 
York, 1994). 

[2] Y.T. Lau, J.M. Finn, and E. Ott, Phys. Rev. Lett. 66 (1991) 978. 

[3] Q. Chen, M. Ding, and E. Ott, Phys. Lett. A 145 (1990) 93. 

[4] K. Sandau, Physica A 233 (1996) 1. 

[5] H.E. Nusse and J.A. Yorke, Physica D 36 (1989) 137. 

[6] V.J. Martinez, R. Dominguez-Tenreiro, and L.J. Roy, Phys. Rev. E 47 
(1993) 735. 

[7] S. Bleher, E. Ott, and C. Grebogi, Phys. Rev. Lett. 63 (1989) 919. 



8 



Figure 1: Cnior R — 0.37: Ci (light grey), C2 (grey), C3 (dark grey) and C4 
(black). 



Figure 2: In Kn for i? = 0.37: n = 5 (least inclined curve), n = 6, n = 7 and 
n = 8 (most inclined curve). 



Figure 3: Estimates of the dimension of fl^^^Cm as a function of R: Dh 
(boxes), Dc (stars) and Dg (triangles). The size of the boxes, stars and 
triangles corresponds approximately to the statistical uncertainty of the di- 
mension. 



Figure 4: Portrait of the basins of the tetrahedron scatterer for: (a) R = 0.48; 
(b) R — 0.37. The initial conditions were chosen on a grid of 400 x 400, on 

the plane P. Regions in black and white correspond to orbits that escape to 
X +00 and X —>■ —00, respectively. Regions in grey correspond to orbits 
not scattered. 
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